Non-canonical functions of UHRF1 maintain DNA methylation homeostasis in cancer cells.

DNA methylation is an essential epigenetic chromatin modification, and its maintenance in mammals requires the protein UHRF1. It is yet unclear if UHRF1 functions solely by stimulating DNA methylation maintenance by DNMT1, or if it has important additional functions. Using degron alleles, we show that UHRF1 depletion causes a much greater loss of DNA methylation than DNMT1 depletion. This is not caused by passive demethylation as UHRF1-depleted cells proliferate more slowly than DNMT1-depleted cells. Instead, bioinformatics, proteomics and genetics experiments establish that UHRF1, besides activating DNMT1, interacts with DNMT3A and DNMT3B and promotes their activity. In addition, we show that UHRF1 antagonizes active DNA demethylation by TET2. Therefore, UHRF1 has non-canonical roles that contribute importantly to DNA methylation homeostasis; these findings have practical implications for epigenetics in health and disease.

RepEnTools: an automated repeat enrichment analysis package for ChIP-seq data reveals hUHRF1 Tandem-Tudor domain enrichment in young repeats.

BACKGROUND: Repeat elements (REs) play important roles for cell function in health and disease. However, RE enrichment analysis in short-read high-throughput sequencing (HTS) data, such as ChIP-seq, is a challenging task. RESULTS: Here, we present RepEnTools, a software package for genome-wide RE enrichment analysis of ChIP-seq and similar chromatin pulldown experiments. Our analysis package bundles together various software with carefully chosen and validated settings to provide a complete solution for RE analysis, starting from raw input files to tabular and graphical outputs. RepEnTools implementations are easily accessible even with minimal IT skills (Galaxy/UNIX). To demonstrate the performance of RepEnTools, we analysed chromatin pulldown data by the human UHRF1 TTD protein domain and discovered enrichment of TTD binding on young primate and hominid specific polymorphic repeats (SVA, L1PA1/L1HS) overlapping known enhancers and decorated with H3K4me1-K9me2/3 modifications. We corroborated these new bioinformatic findings with experimental data by qPCR assays using newly developed primate and hominid specific qPCR assays which complement similar research tools. Finally, we analysed mouse UHRF1 ChIP-seq data with RepEnTools and showed that the endogenous mUHRF1 protein colocalizes with H3K4me1-H3K9me3 on promoters of REs which were silenced by UHRF1. These new data suggest a functional role for UHRF1 in silencing of REs that is mediated by TTD binding to the H3K4me1-K9me3 double mark and conserved in two mammalian species. CONCLUSIONS: RepEnTools improves the previously available programmes for RE enrichment analysis in chromatin pulldown studies by leveraging new tools, enhancing accessibility and adding some key functions. RepEnTools can analyse RE enrichment rapidly, efficiently, and accurately, providing the community with an up-to-date, reliable and accessible tool for this important type of analysis.

Modular dual-color BiAD sensors for locus-specific readout of epigenome modifications in single cells.

Dynamic changes in the epigenome at defined genomic loci play crucial roles during cellular differentiation and disease development. Here, we developed dual-color bimolecular anchor detector (BiAD) sensors for high-sensitivity readout of locus-specific epigenome modifications by fluorescence microscopy. Our BiAD sensors comprise an sgRNA/dCas9 complex as anchor and double chromatin reader domains as detector modules, both fused to complementary parts of a split IFP2.0 fluorophore, enabling its reconstitution upon binding of both parts in close proximity. In addition, a YPet fluorophore is recruited to the sgRNA to mark the genomic locus of interest. With these dual-color BiAD sensors, we detected H3K9me2/3 and DNA methylation and their dynamic changes upon RNAi or inhibitor treatment with high sensitivity at endogenous genomic regions. Furthermore, we showcased locus-specific H3K36me2/3 readout as well as H3K27me3 and H3K9me2/3 enrichment on the inactive X chromosome, highlighting the broad applicability of our dual-color BiAD sensors for single-cell epigenome studies.

Repeat DNA methylation is modulated by adherens junction signaling.

Through its involvement in gene transcription and heterochromatin formation, DNA methylation regulates how cells interact with their environment. Nevertheless, the extracellular signaling cues that modulate the distribution of this central chromatin modification are largely unclear. DNA methylation is highly abundant at repetitive elements, but its investigation in live cells has been complicated by methodological challenges. Utilizing a CRISPR/dCas9 biosensor that reads DNA methylation of human α-satellite repeats in live cells, we here uncover a signaling pathway linking the chromatin and transcriptional state of repetitive elements to epithelial adherens junction integrity. Specifically, we find that in confluent breast epithelial cell monolayers, α-satellite repeat methylation is reduced by comparison to low density cultures. This is coupled with increased transcriptional activity at repeats. Through comprehensive perturbation experiments, we identify the junctional protein E-cadherin, which links to the actin cytoskeleton, as a central molecular player for signal relay into the nucleus. Furthermore, we find that this pathway is impaired in cancer cells that lack E-cadherin and are not contact-inhibited. This suggests that the molecular connection between cell density and repetitive element methylation could play a role in the maintenance of epithelial tissue homeostasis.

Locus-Specific and Stable DNA Demethylation at the H19/IGF2 ICR1 by Epigenome Editing Using a dCas9-SunTag System and the Catalytic Domain of TET1.

DNA methylation is critically involved in the regulation of chromatin states and cell-type-specific gene expression. The exclusive expression of imprinted genes from either the maternal or the paternal allele is regulated by allele-specific DNA methylation at imprinting control regions (ICRs). Aberrant DNA hyper- or hypomethylation at the ICR1 of the H19/IGF2 imprinting locus is characteristic for the imprinting disorders Beckwith-Wiedemann syndrome (BWS) and Silver-Russell syndrome (SRS), respectively. In this paper, we performed epigenome editing to induce targeted DNA demethylation at ICR1 in HEK293 cells using dCas9-SunTag and the catalytic domain of TET1. 5-methylcytosine (5mC) levels at the target locus were reduced up to 90% and, 27 days after transient transfection, >60% demethylation was still observed. Consistent with the stable demethylation of CTCF-binding sites within the ICR1, the occupancy of the DNA methylation-sensitive insulator CTCF protein increased by >2-fold throughout the 27 days. Additionally, the H19 expression was increased by 2-fold stably, while IGF2 was repressed though only transiently. Our data illustrate the ability of epigenome editing to implement long-term changes in DNA methylation at imprinting control regions after a single transient treatment, potentially paving the way for therapeutic epigenome editing approaches in the treatment of imprinting disorders.

Targeting oncogenic TERT promoter variants by allele-specific epigenome editing.

BACKGROUND: Activation of dominant oncogenes by small or structural genomic alterations is a common driver mechanism in many cancers. Silencing of such dominantly activated oncogenic alleles, thus, is a promising strategy to treat cancer. Recently, allele-specific epigenome editing (ASEE) has been described as a means to reduce transcription of genes in an allele-specific manner. In cancer, specificity to an oncogenic allele can be reached by either targeting directly a pathogenic single-nucleotide variant or a polymorphic single-nucleotide variant linked to the oncogenic allele. To investigate the potential of ASEE in cancer, we here explored this approach by targeting variants at the TERT promoter region. The TERT promoter region has been described as one of the most frequently mutated non-coding cancer drivers. RESULTS: Sequencing of the TERT promoter in cancer cell lines showed 53% (41/77) to contain at least one heterozygous sequence variant allowing allele distinction. We chose the hepatoblastoma cell line Hep-G2 and the lung cancer cell line A-549 for this proof-of-principle study, as they contained two different kinds of variants, namely the activating mutation C228T in the TERT core promoter and the common SNP rs2853669 in the THOR region, respectively. These variants were targeted in an allele-specific manner using sgRNA-guided dCas9-DNMT3A-3L complexes. In both cell lines, we successfully introduced DNA methylation specifically to the on-target allele of the TERT promoter with limited background methylation on the off-target allele or an off-target locus (VEGFA), respectively. We observed a maximum CpG methylation gain of 39% and 76% on the target allele when targeting the activating mutation and the common SNP, respectively. The epigenome editing translated into reduced TERT RNA expression in Hep-G2. CONCLUSIONS: We applied an ASEE-mediated approach to silence TERT allele specifically. Our results show that the concept of dominant oncogene inactivation by allele-specific epigenome editing can be successfully translated into cancer models. This new strategy may have important advantages in comparison with existing therapeutic approaches, e.g., targeting telomerase, especially with regard to reducing adverse side effects.

Development of super-specific epigenome editing by targeted allele-specific DNA methylation.

BACKGROUND: Epigenome editing refers to the targeted reprogramming of genomic loci using an EpiEditor which may consist of an sgRNA/dCas9 complex that recruits DNMT3A/3L to the target locus. Methylation of the locus can lead to a modulation of gene expression. Allele-specific DNA methylation (ASM) refers to the targeted methylation delivery only to one allele of a locus. In the context of diseases caused by a dominant mutation, the selective DNA methylation of the mutant allele could be used to repress its expression but retain the functionality of the normal gene. RESULTS: To set up allele-specific targeted DNA methylation, target regions were selected from hypomethylated CGIs bearing a heterozygous SNP in their promoters in the HEK293 cell line. We aimed at delivering maximum DNA methylation with highest allelic specificity in the targeted regions. Placing SNPs in the PAM or seed regions of the sgRNA, we designed 24 different sgRNAs targeting single alleles in 14 different gene loci. We achieved efficient ASM in multiple cases, such as ISG15, MSH6, GPD1L, MRPL52, PDE8A, NARF, DAP3, and GSPT1, which in best cases led to five to tenfold stronger average DNA methylation at the on-target allele and absolute differences in the DNA methylation gain at on- and off-target alleles of > 50%. In general, loci with the allele discriminatory SNP positioned in the PAM region showed higher success rate of ASM and better specificity. Highest DNA methylation was observed on day 3 after transfection followed by a gradual decline. In selected cases, ASM was stable up to 11 days in HEK293 cells and it led up to a 3.6-fold change in allelic expression ratios. CONCLUSIONS: We successfully delivered ASM at multiple genomic loci with high specificity, efficiency and stability. This form of super-specific epigenome editing could find applications in the treatment of diseases caused by dominant mutations, because it allows silencing of the mutant allele without repression of the expression of the normal allele thereby minimizing potential side-effects of the treatment.

Refined read-out: The hUHRF1 Tandem-Tudor domain prefers binding to histone H3 tails containing K4me1 in the context of H3K9me2/3.

UHRF1 is an essential chromatin protein required for DNA methylation maintenance, mammalian development, and gene regulation. We investigated the Tandem-Tudor domain (TTD) of human UHRF1 that is known to bind H3K9me2/3 histones and is a major driver of UHRF1 localization in cells. We verified binding to H3K9me2/3 but unexpectedly discovered stronger binding to H3 peptides and mononucleosomes containing K9me2/3 with additional K4me1. We investigated the combined binding of TTD to H3K4me1-K9me2/3 versus H3K9me2/3 alone, engineered mutants with specific and differential changes of binding, and discovered a novel read-out mechanism for H3K4me1 in an H3K9me2/3 context that is based on the interaction of R207 with the H3K4me1 methyl group and on counting the H-bond capacity of H3K4. Individual TTD mutants showed up to a 10,000-fold preference for the double-modified peptides, suggesting that after a conformational change, WT TTD could exhibit similar effects. The frequent appearance of H3K4me1-K9me2 regions in human chromatin demonstrated in our TTD chromatin pull-down and ChIP-western blot data suggests that it has specific biological roles. Chromatin pull-down of TTD from HepG2 cells and full-length murine UHRF1 ChIP-seq data correlate with H3K4me1 profiles indicating that the H3K4me1-K9me2/3 interaction of TTD influences chromatin binding of full-length UHRF1. We demonstrate the H3K4me1-K9me2/3 specific binding of UHRF1-TTD to enhancers and promoters of cell-type-specific genes at the flanks of cell-type-specific transcription factor binding sites, and provided evidence supporting an H3K4me1-K9me2/3 dependent and TTD mediated downregulation of these genes by UHRF1. All these findings illustrate the important physiological function of UHRF1-TTD binding to H3K4me1-K9me2/3 double marks in a cellular context.

The T1150A cancer mutant of the protein lysine dimethyltransferase NSD2 can introduce H3K36 trimethylation.

Protein lysine methyltransferases (PKMTs) play essential roles in gene expression regulation and cancer development. Somatic mutations in PKMTs are frequently observed in cancer cells. In biochemical experiments, we show here that the NSD1 mutations Y1971C, R2017Q, and R2017L observed mostly in solid cancers are catalytically inactive suggesting that NSD1 acts as a tumor suppressor gene in these tumors. In contrast, the frequently observed T1150A in NSD2 and its T2029A counterpart in NSD1, both observed in leukemia, are hyperactive and introduce up to three methyl groups in H3K36 in biochemical and cellular assays, while wildtype NSD2 and NSD1 only introduce up to two methyl groups. In Molecular Dynamics simulations, we determined key mechanistic and structural features controlling the product specificity of this class of enzymes. Simulations with NSD2 revealed that H3K36me3 formation is possible due to an enlarged active site pocket of T1150A and loss of direct contacts of T1150 to critical residues which regulate the product specificity of NSD2. Bioinformatic analyses of published data suggested that the generation of H3K36me3 by NSD2 T1150A could alter gene regulation by antagonizing H3K27me3 finally leading to the upregulation of oncogenes.

On the accuracy of the epigenetic copy machine: comprehensive specificity analysis of the DNMT1 DNA methyltransferase.

The specificity of DNMT1 for hemimethylated DNA is a central feature for the inheritance of DNA methylation. We investigated this property in competitive methylation kinetics using hemimethylated (HM), hemihydroxymethylated (OH) and unmethylated (UM) substrates with single CpG sites in a randomized sequence context. DNMT1 shows a strong flanking sequence dependent HM/UM specificity of 80-fold on average, which is slightly enhanced on long hemimethylated DNA substrates. To explain this strong effect of a single methyl group, we propose a novel model in which the presence of the 5mC methyl group changes the conformation of the DNMT1-DNA complex into an active conformation by steric repulsion. The HM/OH preference is flanking sequence dependent and on average only 13-fold, indicating that passive DNA demethylation by 5hmC generation is not efficient in many flanking contexts. The CXXC domain of DNMT1 has a moderate flanking sequence dependent contribution to HM/UM specificity during DNA association to DNMT1, but not if DNMT1 methylates long DNA molecules in processive methylation mode. Comparison of genomic methylation patterns from mouse ES cell lines with various deletions of DNMTs and TETs with our data revealed that the UM specificity profile is most related to cellular methylation patterns, indicating that de novo methylation activity of DNMT1 shapes the DNA methylome in these cells.

Efficient Targeted DNA Methylation with dCas9-Coupled DNMT3A-DNMT3L Methyltransferase.

Epigenome editing is a powerful approach for the establishment of a chromatin environment with desired properties at a selected genomic locus, which is used to influence the transcription of target genes and to study properties and functions of gene regulatory elements. Targeted DNA methylation is one of the most often used types of epigenome editing, which typically aims for gene silencing by methylation of gene promoters. Here, we describe the design principles of EpiEditors for targeted DNA methylation and provide step-by-step guidelines for the realization of this approach. We focus on the dCas9 protein as the state-of-the-art DNA targeting module fused to 10×SunTag as the most frequently used system for editing enhancement. Further, we discuss different flavors of DNA methyltransferase modules used for this purpose including the most specific variants developed recently. Finally, we explain the principles of gRNA selection, outline the setup of the cell culture experiments, and briefly introduce the available options for the downstream DNA methylation data analysis.

DNA methyltransferase DNMT3A forms interaction networks with the CpG site and flanking sequence elements for efficient methylation.

Specific DNA methylation at CpG and non-CpG sites is essential for chromatin regulation. The DNA methyltransferase DNMT3A interacts with target sites surrounded by variable DNA sequences with its TRD and RD loops, but the functional necessity of these interactions is unclear. We investigated CpG and non-CpG methylation in a randomized sequence context using WT DNMT3A and several DNMT3A variants containing mutations at DNA-interacting residues. Our data revealed that the flanking sequence of target sites between the -2 and up to the +8 position modulates methylation rates >100-fold. Non-CpG methylation flanking preferences were even stronger and favor C(+1). R836 and N838 in concert mediate recognition of the CpG guanine. R836 changes its conformation in a flanking sequence-dependent manner and either contacts the CpG guanine or the +1/+2 flank, thereby coupling the interaction with both sequence elements. R836 suppresses activity at CNT sites but supports methylation of CAC substrates, the preferred target for non-CpG methylation of DNMT3A in cells. N838 helps to balance this effect and prevent the preference for C(+1) from becoming too strong. Surprisingly, we found L883 reduces DNMT3A activity despite being highly conserved in evolution. However, mutations at L883 disrupt the DNMT3A-specific DNA interactions of the RD loop, leading to altered flanking sequence preferences. Similar effects occur after the R882H mutation in cancer cells. Our data reveal that DNMT3A forms flexible and interdependent interaction networks with the CpG guanine and flanking residues that ensure recognition of the CpG and efficient methylation of the cytosine in contexts of variable flanking sequences.

Methylation of recombinant mononucleosomes by DNMT3A demonstrates efficient linker DNA methylation and a role of H3K36me3.

Recently, the structure of the DNMT3A2/3B3 heterotetramer complex bound to a mononucleosome was reported. Here, we investigate DNA methylation of recombinant unmodified, H3KC4me3 and H3KC36me3 containing mononucleosomes by DNMT3A2, DNMT3A catalytic domain (DNMT3AC) and the DNMT3AC/3B3C complex. We show strong protection of the nucleosomal bound DNA against methylation, but efficient linker-DNA methylation next to the nucleosome core. High and low methylation levels of two specific CpG sites next to the nucleosome core agree well with details of the DNMT3A2/3B3-nucleosome structure. Linker DNA methylation next to the nucleosome is increased in the absence of H3K4me3, likely caused by binding of the H3-tail to the ADD domain leading to relief of autoinhibition. Our data demonstrate a strong stimulatory effect of H3K36me3 on linker DNA methylation, which is independent of the DNMT3A-PWWP domain. This observation reveals a direct functional role of H3K36me3 on the stimulation of DNA methylation, which could be explained by hindering the interaction of the H3-tail and the linker DNA. We propose an evolutionary model in which the direct stimulatory effect of H3K36me3 on DNA methylation preceded its signaling function, which could explain the evolutionary origin of the widely distributed "active gene body-H3K36me3-DNA methylation" connection.

Preferential Self-interaction of DNA Methyltransferase DNMT3A Subunits Containing the R882H Cancer Mutation Leads to Dominant Changes of Flanking Sequence Preferences.

Somatic R882H DNMT3A mutations occur frequently in AML, but their pathogenic mechanism is unclear. As R882H mutations usually are heterozygous, wildtype (WT) and R882H subunits co-exist in affected cells. R882 is located in the RD interface of DNMT3A tetramers, which forms the DNA binding site. R882H causes strong changes in the flanking sequence preferences of DNMT3A. Here, we analyzed flanking sequence preferences for CGNNNN sites showing that most disfavored sites are methylated 4-5 fold slower by R882H than WT, while it methylates most preferred sites 2-fold faster. Overall, R882H was more active than WT at 13% and less active at 52% of all CGNNNN sites. We prepared mixed DNMT3A heterotetramers containing WT and R882H subunits and show that mixed complexes preferentially assemble with an R882H/R882H RD interface. Structural comparisons and MD simulations confirmed the conclusion that the R882H RD interface is more stable than that of WT, in part because H882 forms an inter-subunit contact in the RD interface, while R882 contacts the DNA. As the subunits at the RD interface contribute the two active centers to the DNMT3A tetramer, R882H characteristic flanking sequence preferences of DNMT3A were observed in mixed tetrameric complexes containing WT and R882H subunits, and they are not diluted by the "averaged" effects of mixed or WT interfaces. Hence, R882H has a dominant effect on the flanking sequence preferences and other catalytic properties of DNMT3A in samples containing WT and R882H subunits, which may explain its pathogenic effect in heterozygous state.

Flanking sequences influence the activity of TET1 and TET2 methylcytosine dioxygenases and affect genomic 5hmC patterns.

TET dioxygenases convert 5-methylcytosine (5mC) preferentially in a CpG context into 5-hydroxymethylcytosine (5hmC) and higher oxidized forms, thereby initiating DNA demethylation, but details regarding the effects of the DNA sequences flanking the target 5mC site on TET activity are unknown. We investigated oxidation of libraries of DNA substrates containing one 5mC or 5hmC residue in randomized sequence context using single molecule readout of oxidation activity and sequence and show pronounced 20 and 70-fold flanking sequence effects on the catalytic activities of TET1 and TET2, respectively. Flanking sequence preferences were similar for TET1 and TET2 and also for 5mC and 5hmC substrates. Enhanced flanking sequence preferences were observed at non-CpG sites together with profound effects of flanking sequences on the specificity of TET2. TET flanking sequence preferences are reflected in genome-wide and local patterns of 5hmC and DNA demethylation in human and mouse cells indicating that they influence genomic DNA modification patterns in combination with locus specific targeting of TET enzymes.

Deep Enzymology Studies on DNA Methyltransferases Reveal Novel Connections between Flanking Sequences and Enzyme Activity.

DNA interacting enzymes recognize their target sequences embedded in variable flanking sequence context. The influence of flanking sequences on enzymatic activities of DNA methyltransferases (DNMTs) can be systematically studied with "deep enzymology" approaches using pools of double-stranded DNA substrates, which contain target sites in random flanking sequence context. After incubation with DNMTs and bisulfite conversion, the methylation states and flanking sequences of individual DNA molecules are determined by NGS. Deep enzymology studies with different human and mouse DNMTs revealed strong influences of flanking sequences on their CpG and non-CpG methylation activity and the structures of DNMT-DNA complexes. Differences in flanking sequence preferences of DNMT3A and DNMT3B were shown to be related to the prominent role of DNMT3B in the methylation of human SATII repeat elements. Mutational studies in DNMT3B discovered alternative interaction networks between the enzyme and the DNA leading to a partial equalization of the effects of different flanking sequences. Structural studies in DNMT1 revealed striking correlations between enzymatic activities and flanking sequence dependent conformational changes upon DNA binding. Correlation of the biochemical data with cellular methylation patterns demonstrated that flanking sequence preferences are an important parameter that influences genomic DNA methylation patterns together with other mechanisms targeting DNMTs to genomic sites.

Structural and biochemical insight into the mechanism of dual CpG site binding and methylation by the DNMT3A DNA methyltransferase.

DNMT3A/3L heterotetramers contain two active centers binding CpG sites at 12 bp distance, however their interaction with DNA not containing this feature is unclear. Using randomized substrates, we observed preferential co-methylation of CpG sites with 6, 9 and 12 bp spacing by DNMT3A and DNMT3A/3L. Co-methylation was favored by AT bases between the 12 bp spaced CpG sites consistent with their increased bending flexibility. SFM analyses of DNMT3A/3L complexes bound to CpG sites with 12 bp spacing revealed either single heterotetramers inducing 40° DNA bending as observed in the X-ray structure, or two heterotetramers bound side-by-side to the DNA yielding 80° bending. SFM data of DNMT3A/3L bound to CpG sites spaced by 6 and 9 bp revealed binding of two heterotetramers and 100° DNA bending. Modeling showed that for 6 bp distance between CpG sites, two DNMT3A/3L heterotetramers could bind side-by-side on the DNA similarly as for 12 bp distance, but with each CpG bound by a different heterotetramer. For 9 bp spacing our model invokes a tetramer swap of the bound DNA. These additional DNA interaction modes explain how DNMT3A and DNMT3A/3L overcome their structural preference for CpG sites with 12 bp spacing during the methylation of natural DNA.

Genome-wide investigation of the dynamic changes of epigenome modifications after global DNA methylation editing.

Chromatin properties are regulated by complex networks of epigenome modifications. Currently, it is unclear how these modifications interact and if they control downstream effects such as gene expression. We employed promiscuous chromatin binding of a zinc finger fused catalytic domain of DNMT3A to introduce DNA methylation in HEK293 cells at many CpG islands (CGIs) and systematically investigated the dynamics of the introduced DNA methylation and the consequent changes of the epigenome network. We observed efficient methylation at thousands of CGIs, but it was unstable at about 90% of them, highlighting the power of genome-wide molecular processes that protect CGIs against DNA methylation. Partially stable methylation was observed at about 1000 CGIs, which showed enrichment in H3K27me3. Globally, the introduced DNA methylation strongly correlated with a decrease in gene expression indicating a direct effect. Similarly, global but transient reductions in H3K4me3 and H3K27ac were observed after DNA methylation but no changes were found for H3K9me3 and H3K36me3. Our data provide a global and time-resolved view on the network of epigenome modifications, their connections with DNA methylation and the responses triggered by artificial DNA methylation revealing a direct repressive effect of DNA methylation in CGIs on H3K4me3, histone acetylation, and gene expression.

Complex DNA sequence readout mechanisms of the DNMT3B DNA methyltransferase.

DNA methyltransferases interact with their CpG target sites in the context of variable flanking sequences. We investigated DNA methylation by the human DNMT3B catalytic domain using substrate pools containing CpX target sites in randomized flanking context and identified combined effects of CpG recognition and flanking sequence interaction together with complex contact networks involved in balancing the interaction with different flanking sites. DNA methylation rates were more affected by flanking sequences at non-CpG than at CpG sites. We show that T775 has an essential dynamic role in the catalytic mechanism of DNMT3B. Moreover, we identify six amino acid residues in the DNA-binding interface of DNMT3B (N652, N656, N658, K777, N779, and R823), which are involved in the equalization of methylation rates of CpG sites in favored and disfavored sequence contexts by forming compensatory interactions to the flanking residues including a CpG specific contact to an A at the +1 flanking site. Non-CpG flanking preferences of DNMT3B are highly correlated with non-CpG methylation patterns in human cells. Comparison of the flanking sequence preferences of human and mouse DNMT3B revealed subtle differences suggesting a co-evolution of flanking sequence preferences and cellular DNMT targets.